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Abstract. In this work we consider the generahzed Navier-Stoke equations 
with the presence of a damping term in the momentum equation. The prob- 
lem studied here derives from the set of equations which govern the isothermal 
flow of incompressible, homogeneous and non-Newtonian fluids. For the gen- 
eralized Navier-Stokes problem with damping, we prove the existence of weak 
solutions by using regularization techniques, the theory of monotone operators 
and compactness arguments together with the local decomposition of the pres- 
sure and the Lipschitz-truncation method. The existence result proved here 
holds for any q > ^^'^ cr > 1, where q is the exponent of the diffusion 

term and a is the exponent which characterizes the damping term. 
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1. Introduction 

In this work, we shall study the existence of weak solutions for the generalized 
Navier-Stokes equations with damping: 

(1.1) divu = in Qt, 

(1.2) -^-div(|Vu|«"2yu-u®u) +a|u|'""2u^f _yp jj^ Q^. 

supplemented with the following initial and boundary conditions: 

(1.3) u = Uo in O for t = 0, 

(1.4) u = on Tt. 
Here Qt is a general cylinder defined by 

QT:=nx (0, T), with Ft := dn x (0, T), 
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where n C R'^ , N > 2, is a bounded domain with a compact boundary dfl, and 
< T < oo. 

In the scope of Mathematical Fluid Mechanics, u is the velocity field, p stands 
for the pressure divided by the constant density, f is the given forcing term and 
q > 1 is the constant exponent which characterizes the flow. The constant a is 
non-negative and cr > 1 is another constant. 

The damping term q;|u|'^^^u, or sometimes called absorption term, has no direct 
physical justification in Fluid Mechanics, although it might be considered has being 
part of the external body forces field (see There is also a precise theory of 

the absorption of forced plane infinitesimal waves according to the Navier-Stokes 
equations (see |23]). The consideration of damping terms in the generalized Navier- 
Stokes equations it is also useful as a regularization procedure to prove the existence 
of weak solutions for the stationary problems (see [Il|-[14]). At last, but not in last, 
there is also the purely mathematical motivation which goes back to a work about a 
stationary like problem (see [5]), where the authors where mainly interested with the 
important question about compact supported solutions for that problem. During 
the last years, many authors have worked on these kind of modified Navier-Stokes 
type problems, establishing the existence of weak solutions and proving many other 
properties has the uniqueness of weak solutions, their regularity and studying its 
asymptotic behavior. In [lH| we proved the weak solutions of (|l.ip - (|1.4p extinct in 
a finite time for g > 2, provided I < a < 2. This property is well known for the 
generalized Navier-Stokes problem (|l.ip - ()1.4p with a = in the case 1 < g < 2. 
But for (7 > 2 the best one can gets are some decays of fractional and exponential 
order (see e.g. [7]). In [4] we have studied the problem (|l.l[) - (|1.4p in the particular 
case of g = 2. There, we have proved the existence of weak solutions, its uniqueness 
and some asymptotic properties. We carried out an analogous study in j5] for the 
Oberbeck-Boussinesq version of this problem, where besides the usual coupling in 
the buoyancy force, we have considered an extra coupling in the damping term by 
considering a temperature-depending function cr. In |10| the authors have proved 
the existence of weak and strong solutions for the Cauchy problem (|l.ip - (|1.4p in 
R'^ and with q = 2. The damping term is being considered in the context of many 
other physical systems which go from the Shrodinger equations (see e.g. [11]) to 
the Euler equations (see e.g. [19]) and passing by the wave equation (see e.g. |26|). 

With respect to the existence of weak solutions for the original generalized 
Navier-Stokes problem, i.e. (|l.ip - (|1.4p with a = 0, the problem was solved in 
its all full possible (it is open only the case 1 < g < ■^^) extension recently in the 
work [12 • The first existence result to this problem was achieved in [T^] and [T7] for 
q > . Only more or less 40 years later it was possible to improve the existence 
result for lower values of g. In |25| . under the same assumptions of [16] and [17], it 
was proved the existence of weak solutions to the problem (|l.ip - (|1.4p with a = 

for g > max A ht bit earlier to the work [25|, it was proved 

in [23] the existence of a weak solution to the same problem for g > 2^ii. Finally 
in [12] the authors have extended the result [53] to the case g > . It is an open 
problem to prove the existence of weak solutions to the problem (|l.ip - (jl.4p with 
a = Oifl<g< in the case of iV > 2. On the other hand, it seems to be very 
difficult to go bellow the limit g = (for > 2), due to the need of using the 
compact imbedding W^''?(ri) ^ L^(ri). 
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The plan of this work is the following. In Section [T] we introduce the problem 
we shall study here and review some results related with our work. Section [2] is 
devoted to introduce the notation we use throughout the work and to define the 
notion of weak solution we shall consider. Here, we also shall state the main result 
of this paper: Theorem 12.11 The proof of this result is carried out from Section [3] 
to Section [m In Section [T2| we make some remarks about our work, in special, its 
extensions and limitations. 



The notation used throughout this article is largely standard in Mathematical 
Fluid Mechanics - see e.g. [17]. We distinguish tensors and vectors from scalars by 
using boldface letters. For functions and function spaces we will use this distinction 
as well. The symbol C will denote a generic constant - generally a positive one, 
whose value will not be specified; it can change from one inequality to another. 
The dependence of C on other constants or parameters will always be clear from 
the exposition. In this paper, the notations $7 or w stand always for a domain, i.e., 
a connected open subset of R^, iV > 1. By Cn{,oj) we denote the A^-dimensional 
Lesbesgue measure of w. Given fc £ H, we denote by 0*^(0) the space of all k- 
differentiable functions in il. By C^(J7) or X'(ri), we denote the space of all infinity- 
differentiable functions with compact support in fl. The space of distributions 
over 2?(17) is denoted by I?'(r2). If X is a generic Banach space, its dual space is 
denoted by X'. Let 1 < q < oo and il C E,^, with N > 1, be a domain. We 
will use the classical Lebesgue spaces L*(il), whose norm is denoted by || • ||L9(n)- 
For any nonnegative k, W'^'''(J7) denotes the Sobolev space of all functions u E 
L'(n) such that the weak derivatives D"u exist, in the generalized sense, and are 
in L'(ri) for any multi-index a such that < |a| < k. The norm in W''''^{^1) 
is denoted by || • |lw'' «(f2)- The corresponding spaces of vector-valued or tensor- 
valued functions are denoted by boldface letters. All these spaces are Banach 
spaces and the Hilbert framework corresponds to g = 2. In the last case, we use 
the abbreviation W*^'^ = H*^. Given T > and a Banach space X, L'(0, T; X) and 
W'=''?(0,T;X) denote the usual Bochner spaces used in evolutive problems, with 
norms denoted by || • \\hi(o,T;X) and || • || w*'9(o,T;X)- By Cw([0,T];X) we denote the 
subspace of L°°(0, T; X) consisting of functions which are weakly continuous from 
[0,r] into X. 

A very important property satisfied by the tensor S = |Vu|''^^Vu and by the 
damping term |u|°'~^u are expressed in the following lemma which proof we address 
the reader to [6]. 

Lemma 2.1. For all s e (l,oo) and 6 > 0, there exist constants Ci and C2, 
depending on s and N , such that for all ^, 77 € R''^, A'^ > 1, 



In order to define the notion of weak solutions we shall look for, let us introduce 
the usual functional setting of Mathematical Fluid Mechanics: 



2. Weak formulation 



(2.5) 



m-'^ - ivr'v) ■ (e - v) >C2\^- v\'^' m + \v\) 



S-2-S 



(2.6) 



V := {v G C^{n) : divv = 0} ; 



(2.7) 



H := closure of V in 



(2.8) 



:= closure of V in W^^^il) 
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The weak solutions we are interested in are usually called in the sense of Leray-Hopf . 

Definition 2.1. Let N >2 and 1 < q, <t < oo. Assume that Uo G H, f G L^(Qt) 
and let conditions (A)-(D) he fulfilled for any q > \. A vector field \i is a weak 
solution to the problem hl.l\l - ^T7^ , if: 

(1) ue L-'(0,^;H)nL9(0,T;V,)nL-(QT); 
(2) For every (p G C°°{Qt), with div (y9 = and suppi/j CC x [0,T), 

- / u-iftdxdt+ / (|Vu|'^"^Vu- u(g) u) ■.Vtpdxdt + a / \u\"'^'^u ■ (p dxdt 
JQt JQt Jqt 

= / f-ipdxdt+ / uo • v?(0) rfx. 

Jqt JiT. 

Remark 2.1. Note that ^ L'^(fi) for o < q* , where q* is the Sobolev conjugate 
of q, i.e. q* ~ if 1 < q < N , or q* = oo if q > N. Therefore, in this case, we 

look for weak solutions in the class L°°(0, T; H) H i'(0, T; Vq). 

The main result of this work is the following, where it is established the existence 
of weak solutions to the problem (jl.ip - (jl.4p under the minor possible assumptions 
on q and a. We left open only the case of 1 < g < -jf^, for N > 2, which will 
certainly require a different approach. 

Theorem 2.1. Let il be a bounded domain in M,^ , N > 2, with a Lipschitz- 
continuous boundary dfl. Assume that 

(2.9) f gL«'(0,T;V;), 

(2.10) uo G H. 

Then, if 

2N 

(2.11) .>^, 

there exists a weak solution to the problem HLl\) - [l~^\ l, in the sense of Definition \2. 1[ 
for any a > 1. Moreover, any weak solution u G Cw([0,T];H). 

In order to simplify the exposition, we shall assume throughout the rest of this 
work the following simplified assumption of ()2.9p 



(2.12) f = -divF, F G (Qt)- 

The main ingredients of the proof of Theorem l2.1l are the results of the local decom- 
position of the pressure established in |24j and the Lipschitz-truncation method in 
the spirit of [12]. With could also have considered the L°° -truncation method used 
in |24| , but by this method we cannot achieve an existence result for so lower values 
of q as we can with the Lipschitz-truncation method. The proof of Theorem 12.11 
will be the aim of the next sections. 

3. The regularized problem 

We start the proof of Theorem l2.1l bv considering a regularization of the problem 
(|l.ip - p.4p which basically gets rid off the difficulties coming from the convective 
term u(g) u. Let $ G C°°([0, oo)) be a non-increasing function such that < $ < 1 
in [0, cx)), $ = 1 in [0, 1], $ = in [2, oo) and < < 2. For e > 0, we set 

(3.13) $e(s) := $(es), sG[0,oo), 
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(3.18) 



and let us consider the following regularized problem: 
(3.14) divu, = in Qt, 

(3.15) 

-^~dW{\Wu,\''~^\7u,-u,<»u,<^>,{\u,\))+a\u,r-^u, = i-^p,, in Qt, 

(3.16) Ug ~ Uq in O for t = 0, 

(3.17) Ue = on Tt- 

A vector function e L°°(0, T; H) n L'?(0, T; V,) n L^iQr) is a weak solution to 
the problem ([XTi l - fXTT)) . if 

-/ Ue-iptdxdt+ (iVu^l'-^VUe - Ue «) Ue$e(|Ue|)) : V</Jdxdt 

+ a / |Ue|'^"^Ue • (ysdxdt = / F:V</jdxcft+ / Uq • </j(0) fix 

v/Qt Jqt JiT. 

for all <(9 e C°°((3t), with div(^ = and supp(/3 CC Vl x [0,T). 

Proposition 3.1. Let the assumptions of Theorem \2.1\ he fulfilled. Then, for each 
e > 0, there exists a weak solution e L9(0,T; Vg) n Cw([0,r];H) n L'^(Qt) to 
the problem \3.1Ji}( l- ^3.1'T\) . In addition, every weak solution satisfies to the following 
energy equality: 

(3.19) i||u,(t)||^+ / |Vu,|9dx<it+a / |u,|'"dx(it = i||uo||H+ / F : Vu.dxdt 
/or aZHe (0,T). 

Proof. The proof of Proposition 13. II is adapted from the proof of [241 Theorem 3.1]. 
The difference here is the presence of an extra term which results from the damping 
and the aspect of the diffusion term. We shall split this proof into three steps. 

First Step. Let G (0,T] be arbitrarily chosen and let us set 

Mt, := {w e L2(0,T,;H) : ||cp!|l2(o,t.;H) < !}• 

Observing that by the property (|2.5p , the diffusion term is monotonous as well the 
damping term, we can use the theory of monotone operators (c/. |17| Section 2.1] 
together with O Section 9]) to prove that for each w € Mr,, there exists a weak 
solution V G L°°(0,r*;H) n L«(0,r*; V,) n L'"(QtJ to the following system: 

(3.20) divu^O in Qt,, 



(3.21) -^-div(|Vu|9~2Vu)+a|w|'""2i; = f-Vp-div(cc7«)tu$,(|ro|)) in Qt„ 

(3.22) w = Uo in Q. for t = 0, 

(3.23) v = on Vt,. 

Moreover, once the diffusion and damping terms satisfy the monotonicity property 
(|2.5p . the weak solution of p.20p - p.23p is unique. 

Second Step. As a consequence of the previous step, we can define a mapping 

(3.24) K : AfT. L2(0,T,;H) 
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such that to each vj £ Mt, associates the unique weak solution v E L°°(0, T*; H) n 
L''(0,T*;Vg) n L^iQr,)- Testing formaUy ((3?2T|l by the unique weak solution 
V :— K(ti7), with vj € Mt,, integrating over Qt-, with < t < T*, using Young's 
inequality and, at last, the definition of $(:(|tu|), we achieve to 



(3.25) 



II^^IIl-(ot.-h) +Ci / \Vv\''dxdt + C2 [ \v\''dxdt< 

71 +72|kllL2(o,T.;H): 71 := I!uo||h + C3 / |F|«'dxds, 72 := C4 

J Qt.. 



Then setting := min {1/(71 +72),r}, we can prove, from p.25p and due to the 
fact that w G Mt, , that 

(3.26) \m^)\\hio.T,;n) ^ r,(7i +72) < 1 

for all zu e Mt, ■ This proves that K maps Mt, into itself. 

On the other hand, in order to prove the compactness of K, we obtain from 
(IX^ that 

(3.27) ||K(tn)||L,(o,T.;V,) < C^'hi + 72) 

for all w € Mt. ■ Owing to the assumptions (|2.10p and (|2.12p , the right hand side 
of (|3.27p is finite. Then, for the distributive time derivative v' (K(ro))', with 
Tu e Mt, , we can prove that 

lldiv (iV^r-^Vu - tn w<S>,i\w\) - F) ||l,'(o,t.;Vj.) + II h^r"'^^llL-'(Q^j < ^ 

and consequently 

(3.28) w' e Li'(0,T,;V') +L'^'(Q 



T. 



In fact, by virtue of p.25p . it follows the uniform boundedness of |Vu|'' ^Vu in 
L'''(QtJ and of Ivl^-^v in L"' (Qt, )- By assumption ^J^, F € L«'(QrJ- On the 
other hand, using the definition of $e, we can prove that ||tn(X)i:A7$e(|-n7|)|j^^, ^ < 

C'll^llL(o,T.;H)- By ([3:271 and ^M), and once that V, CC H C V;, for q > 
we can apply Aubin-Lions compactness lemma (cf. [21]) to prove that K(Mt, ) is 
relatively compact in L'(0,T*;H). Then, since v G L°°(0,T*;H), by parabolic 
interpolation it follows that ^.{Mt,) is relatively compact in L^(0,T*;H). 
To prove the continuity of K, we consider a sequence Wm in Mt, such that 

— )■ cij in L^(0,T*;H) as m 00. 

By the relative compactness of ^.{Mt, ) in L^(0, T*; H), there exists a subsequence 
vumk such that 

(3.29) K(ro™J->w in L2(0,T,;H), as k ^ 00. 

From the definition of K, the functions Vrrn. K(n7„is.) satisfy to 

- / Vm^ ■ (fit dx.dt + / |Vw,„j^|''^^Vw„ij^ : Vi^dxdt 
Jqt, Jot. 



(3.30) +a / \v,r,J'' ■'vrn^-ipdxdt 

JQt, 



/ (F + tumj^ (g) Ti7m^$e(|w,„J)) : V<^(ix(ii + / Uq • (p(0) dx 



EXISTENCE FOR THE GENERALIZED NAVIER-STOKES EQUATIONS WITH DAMPING 7 



for all ^p E C°°{Qt,), with div (p = and supp<y9 CC x [0,r*). Passing to the 
limit in p.30p by using the appropriated convergence results (see |241 p. 119] and 
[HI p. 236]) and the usual Minty trick (see e.g. [T71 pp. 212-214]), we can prove 
that V = K(tz7). The only difference here is that \vmk\'^~'^'Vnik ^ ^ weakly in 
L'^ (Qt,)-, as fc — > oo. Since Vjru. — v weakly in L'^((3t,), as fc — > oo, there exists 
a subsequence, still denoted by Vm^. , such that — >■ v a.e. in Qt, . In addition, 
because \vmk\'^~'^'^mk is uniformly bounded in L'^ (Qt,), we can apply Lesbesgue's 
theorem of dominated convergence to prove that 

(3.31) IvmX'^Vrn, ^ \v\^-\! strougiy in L"'(QtJ 

and, as a consequence, v = From (|3.29p . we conclude that J<i{zJm) ^ 

K(n7) in L'^(0, T*; H) as m — > oo, which proves the continuity of K. 

Now, applying Schauder's fixed point theorem, there exists a function vt, G Mt, 
such that K(fT, ) = vt, and which is a weak solution to the problem p.l4p - p.l7p 
in the cylinder Qt, . 

Third Step. Testing (|3.15p by the weak solution vt. , integrating over Qt, , proceed- 
ing we did as for p.25p and observing that due to the definition of the term 
resulting from convection is zero, we obtain 

(3.32) IK^T.(n)||^<c(|luo||^ + !lF|l^',,^^^^ 

The estimate p.32p is independent of T* and therefore we can extend vt, as a weak 
solution to the problem p.l4p - ()3.17p in the whole cylinder Qt- 

Finally, the energy relation (|3.19p follows by testing (|3.15p by a weak solution 
and integrating over Qt with < t < T^,. □ 

4. Existence of approximate solutions 

Let u, e L^(0,T; Vq)nL°°(0,T;H)nL'"(QT) be a weak solution to the problem 
([XTi)l - (|XT7)l . From Proposition [O (see ([XTgi) ). we can prove that 

(4.33) I|u.IIl-(ot-h)+ / |Vu,|«dxd<-K /" \u,\^d^dt < C, 

Jqt Jqt 

where, by the assumptions (|2.10p and (|2.12p . C is a positive constant which does 
not depend on e. From (|4.33p we obtain 

(4-34) l|u.|li-(o,T;H) + l|Ue|lL(o,T;V,) < 

(4.35) ||u,||l.(q,) <C. 
Using and it follows that 

(4.36) |||Vu,r-2Vue|lL,'(Q,) <C, 

(4.37) |||u,r-2u,||L.'(Q,) <C. 

On the other hand, by using (|4.34p and the Sobolev imbedding L'^{0, T;Vq) n 
L°°(0,T;H) ^ L«^(Qt) (see [T71 p. 213]), we can prove that 

(4.38) ||u,|| N+2 <C. 

^ ' " "l"— (Qt) - 

As a consequence of ()4.38p and of the definition of $e (see p.l3p ). 

(4.39) ||u,®U,$,(|Ue|)||^,iV+,^^^^ <C. 
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Note that the constants in ()4.34p - ()4.39p are distinct and do not depend on e. From 
()4.34p - (j4.39p . there exists a sequence of positive numbers em such that — > 0, as 

m — > cx), and 

(4.40) u weakly in L'5(0, T; Vq), as m -> oo, 

(4.41) u^^ u weakly in L°"(Qt), as to — j- oo, 

(4.42) iVue^l'^^Vue^^ ^ S weakly in L«'(gT), as m ^ oo, 

(4.43) |Ue^|'^~^Uj^ — > u weakly in L"^' (Qt), as to — > oo, 

N + 2 

(4.44) Ue^ u weakly in L'i~ (0, T; Vq), as m oo, 

(4.45) Ug^$g^ (|ug^ I) — ^ G weakly in L*^ 2iv (Qj'), as to — )■ oo. 

Here we observe that using (|4.43p and arguing as in the proof of Proposition 13.11 
(see p.3ip ). we can prove that 

(4.46) |u,„J'"-2u^^ _^ 1^1^-2^ strongly in L'"'(Qt)- 

As a consequence, we ca write u = |u|'^^^u. Then, using the convergence results 
(|4.40p - (|4.46p . we can pass to the limit £,„ — > in p.lSp with replaced by u^^, 
to obtain 
(4.47) 

-/ u-iftdxdt+a lul" ^^u-ipdxdt+ (S-G-F) : V(y9 dxdt = / Uo-(/3(0)dx 
Jqt Jqt Jqt 

for all ip e C°°{Qt), with div = and supp (p CC x [0, T). 

5. Convergence of the approximated convective term 

In this section we shall prove that G = u. We start by observing that, from 
(|XTS)1 . it follows 

Ug,^ • ift dxdt + a |ue^ |°'~^U(;^ • ipdxdt+ 



(5.48) 

/ {\Vu^J'^-'^Vue^ -u^^(E)u^,^^^^{\u,J) -F) ■.Vipdxdt = 

Jqt 

for aU (f e C^iQr) with div^j = 0. Then, from (j^l^ . and ([T^ . we have 

(5.49) Qg„ := I Vu,„ r-'Vu,,„ - ® Ue„$,„, (|u,„ I) - F e L''(Qt) 
for any r satisfying to 

(5.50) l<.<min{^i^,,' 

Using (|4.37p and (|5.49p - (|5.50p . we can obtain, from (|5.48p . that the distributive 
time derivatives 

(5.51) <^ e L'-(0,T;W-i^^'^(f])) +L'^'(gT). 

Due to the admissible range for r (see ()5.50p ). there always exists a 7 > 1 such that 
the following compact and Sobolev imbeddings hold 

(5.52) Wi;^'(17) CC L^(17) C W-i-(l]), ((r')*)'<7<'Z*, 
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where q* is the Sobolev conjugate of q and r' is the Holder conjugate of r. Then, 
using Aubin-Lions compactness lemma (cf. Simon [21]), we obtain from (|4.40p 
together with (|5.5ip and (|5.52p . and passing to a subsequence, that 

(5.53) -> u strongly in U {0 ,T {Vl)) , as m — > oo. 
Using parabolic interpolation, we obtain from (|4.34p and (|5.53p that 

(5.54) ^ u strongly in L''(0, T; L'*'(f2)) V s : 1 < s < cxj, as m^oo. 

Now, observing that q > is equivalent to q* > 9-^^^, we can choose 7 such 

that q^^ < J < q* such that, in view of (|5.54p . 

N + 2 

u strongly in L''(Qt), s = q — — — , as m — >■ 00. 

In consequence 

(5.55) (g) Ug^^e^^ (juj^ I) — ^ u (X) u strongly in L'tw^ (Qt-), as to —t^ 00. 
Finally, from (|4.45p and ()5.55p , we conclude that G = u (g) u. 

6. Weak continuity 
We start this section by proving that 

(6.56) ue C^([0,r];H). 

We observe that, from (|4.47p . the distributive time derivative Ut is uniquely defined 

by 

(6.57) (ut,^) ^ {dW{S-G)-a\u\^~^u + F,^) V G Cg°(0, T; Y), 
where 

(6.58) Y := V, nV«nL'"(QT)nH, k 
Then we can prove that 

(6.59) Ut G L— (0,r;Y'). 

In fact, due to and immediately follows that divS e L^(0,T; V^) 

and |u|'^~2u ^ L"' (Qt), respectively. By assumption (HHH), F e LT^(0,r;V^). 

That divG e L^{0,T;V'^) follows by Km if fc' ^ q^, which in fact is true 
by our choice of k ( cf. (|6.58p ) . 

Next, let tg G [0,r] be fixed and let tk be a sequence in [0,r] such that 

tk ^ to, as fc — > 00, and such that u{tk) G H V fc G M. 

Then we consider the continuous representant of u in C(0,r;Y'), which exists 
by virtue of (|6.57p and (|6.59p . Finally by means of reflexivity in H and of the 
continuous and dense imbedding of H into Y', we can prove that 

u(tfc) — > u(io) weakly in H, as k 00, 

and whence ()6.56p . 

Now, let us prove that for every t E [0,T] 

(6.60) Ue^(i)— >u(t) weakly in H, as m — > 00. 



q{N + 2)y 



2N 
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Due to (|4.34p . there exists a subsequence u^,^^ (t) e H such that 

Uf^^^ (t) — > T] weakly in H, as m — > oo. 
Arguing as we did for (|6.57p and (|6.59p . the distributive time derivative u^^ e 

N + 1 

L~iv- (0, T; Y') and is uniquely defined by 
(6.61) 

(div(|Vu,„^ r'^n,^^ - u,„^ <g> u,„^ <^,^^ (|u,„ J)) - a|u,„.^ T^'u.^, + F, V') 

for aU (p e Cjf (0,T;Y). In particular, there holds u^^^^ e C(0,r;Y'). Next, we 
introduce rj in (j6.6ip . we use integration by parts and we carry out the passage to 
the limit in the resulting equation by using the convergence results (|4.40p - (|4.45p . 
Combining this equation with the one which results from inserting r/ into (|6.57p and 
integrating by parts, we obtain u{t) = rj, which yields (|6.60p . Finally, combining 
(|6.56p and (|6.60p . we see that also u^^ satisfies to (|6.56p . 

7. Auxiliary results for decomposing the pressure 

Here we make a break in the proof of Theorem 12.11 to show that the results of 
Wolf [12] concerned with the local decomposition of the pressure still hold in the 
case of the momentum equation modified by the presence of the damping term. 
For, let uj' be a fixed but arbitrary open bounded subset of 51 such that 

(7.62) uj' CC n and duj' e 

Given s such that 1 < s < cxj, lets us consider the following auxiliary function 
spaces related with the Helmholtz-Weyl decomposition (c/. [24l Section 2], see also 
[111 Section III.l] and [ini): 

(7.63) A-^(w') {a e L^(w') : a = A/, / e wI''{lu')}; 

{7M)B%uj') := l^beB'iuj') : J 6dx = o|, B'' {uj') := {b e {w') : Ab ^ 0} . 

Proposition 7.1. Let Q e L''i(cj^), q E h"^ (w'rp) , with 1 < si, S2 < oo, and 

(7.65) uG Cw([0,T];H) 
where here H is defined over lj' . Suppose that 

(7.66) — / u ■ ift dxdt + Q : Vip dxdt + / q_- f dxdt = 

for all ip € Cg°(a;^) with diviy9 — and where lo'j. = oj' x (0,T) and ui' satisfies to 
i7.6^) . Then there exist unique functions 

(7.67) / e L'*°(0,T;A"«(a;')), 

(7.68) / e Cw([0,r];B-^"(^')), 
where sq can be taken such that 

(7.69) 1 < So < inin{si, S2} , 
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such that 



(7.70) 



u • ipt d^dt + 
p^^Abfip dxdt 



Q : V(/3 dxdt 



P 



dt 



dxdt 



q • tp dxdt 

u(0) • ip{Q)dt 



for all ip G C°°(w^), withsuppip CC w'x[0,r). In addition, the following estimates 
hold 



(7.71) 



\\P' 



|L'-o(a;^) < Ci (^||Q||L=i(tui,) + ||q||L--2(cu^)) , 



(7.72) ||p''||l~(0,T;L=o(c^')) < C2 (||u||l-(0,T;L2(c^')) + IIQIlL = i(u;i,) + l|q||L=2(c^; 

where Ci and C2 are positive constants depending only on Si (i = 0, l,2j, N and 

Proof. Let t/j G Cg°(a;') with div?/; and let 77 e Cg°(0, T). Inserting ip = tprj into 
(|7.66p and using Fubini's theorem, we obtain 



Prjdt + / 7?7 dt, 



a?7' dt - 

Jo 

where for t e [0,T] 

a{t):= / u(i)-?Adx, /?(*):= / Q(t):VV'rfx, 7(t) := / q(t)-V'dx. 

J Uj' J Uj' J U)' 

Since Q e L"i(w^) and q £ L^2(w^), we have {3 £ L'*i(0,r) and 7 e L^^^q^T). 
In consequence, a G W^'''''(a;^) for any sq such that 1 < Sq < min{si,S2}- By 
Sobolev's imbedding theorem, a is represented by a continuous function, which we 
still denote by a. Using integration by parts, we can represent 



a{t)^a{0)+ I3{s)ds+ 7(5) ds Vte(0,r) 



(7.73) 

Let t € (0,r) be arbitrarily chosen. Using Fubini's theorem, the identity (|7.73[) 
reads 



iu{t) - u(0)) • i' + q(t) • ij + Q(t) : VV' 



dx = 0, 



Q(t) := / Q(s)ds, q(t) / q(.s) 



where 



Now, by the results of M.E. Bogovskii (c/. [151 Theorem III.3.1]) and of K.I. Pilet- 
skas (c/. p31 Theorem III. 5. 2]), there exists a unique function 



p(t) G L"»(w') with / p{t)dx^Q 



such that 
(7.74) 



(u(t) - u(0)) • V + q(t) • V' + Q(i) : 



dx 



p{t)dYvip dx 



for all Ip e Wj''*"(w'). In addition, 
(7.75) \\p\t)Us,^^,^<c{\\n{t)-u{ 



1l=o(<^') + !|Q(t)|lL-'o(a;') + ||q(t)||L=o(c^'), 
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On the other hand, by the apphcation of Hehnhohz-Weyl decomposition of L'*" (w') 
(c/. [201 Theorem 1.4]), there exist p°{t) e A'^'iuj') and e B^'-'iaj') such that 

(7.76) p(0=/W+/W in ^' 

and where A''(aj') and B''(w') are defined in (|7.63p - (|7.64p . Moreover, the sum 
A^" {uj') + B''° (w') is direct. Now, from (|7.74[) and, as a consequence of the assump- 
tion (|7.65p . we can infer that 

(7.77) peC„([0,r];L^°(w')). 
From (|776l) and ((7J7| it follows that 

(7.78) p"eCw([0,r];A^«(c^'))- 

(7.79) / e Cw([0,r];B-^»(c^'))- 

As a consequence of (|7.75l) . (|7.76p . we can derive (j7.72p . Moreover, inserting -0 = V(/) 
in (|7.74p . for G C3°(a;'), using (|7.76p . integrating by parts the resulting equation 
and observing that, in view of (|7.79p . Ap'' — 0, and, in of view of (|7.65p . divu = 0, 
we obtain 

(7.80) / (q(<) • V(/) + Q(i) ; V^^) dx = /" f{t)A(j)dx 

for all (f) G Cjf (oj'). Now. using (|7.80p and proceeding as in [211 pp. 115-116], we 
can prove that / G W^'"" (0, T; A"°(w')) and 



~dt 



(7.81) ||/||l=o(^^) <c(||Q||L=i(„^) + ||q||L-2(^^)), where 

and the constant C depends only on Sj {i = 1, 2, 3), iV and w^. Whence (|7.7ip holds. 
Finally, the identity (|7.70p follows by integrating (|7.74p over (0, T), replacing there 
(/? by If and using ((7J6)) . (fTTSl . (TTTQI and the definition of p° given in (fTST]) . 
The uniqueness of p*^ and p'* follow from ()7.7ip and (j7.72p , respectively. □ 

8. Decomposition of the pressure. 

Let us continue with the proof of Theorem 12. II Using the results of the previous 
section, we shall decompose the pressure into two different components. For, let lo' 
be a domain in the conditions of the previous section (see ()7.62p ). Clearly, in view 
of (|5.48p and with the notation introduced in (|5.49p . we can write 

(8.82) -/ u^^ ■ ft dndt + a \u^J''-'^u^^-Lpdxdt+ Q^^^ : ip dxdt = 

for all ip e C^Wt) with div = and where uj'j, := uj' x (0, T). The resuhs (18^82)) 

and ()6.56p allow us to apply Proposition 17.11 with Q = Qe„^, q = ajue,^ |'^~^Ue^, 
si = r, S2 ~ o' and 

(8.83) So = ro := min{r, ct'}. 

Observe that by (|5.50p . r < 2 and consequently < 2. Therefore we can say that 
exist unique functions 

(8.84) pO^^ e L''°(0,r;A'^"(w')), 

(8.85) eCw([0,r];B'^°(w')), 
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such that 

— / Ue ■iptdx.dt + a I |uc I'^^'^Ue ■ipdiidt+ I Qe ■.Vipdxdt = 

(8.86) /"^ " \ " 
/ p° dwip dxdt - / p'^' -^-^ dxdt + Uq ■ <^(0) dt 

for aU If £ C°°{u!'rp), with supp^s CC uj' x [0,r). In addition, by the same resuh, 
the following estimates hold 

(8.87) lbLllL'-o(-^) < (||Q.„J|l^(.^) + i|u.,J|L.'(„^)) , 



lL°°(0,T;L'-o(cj')) 



< C2 (l|Ue,J|L~(0,T;H) + IIQe„llL'-(i^^) + 1 1 Ue„J | l-' (c^^ 



where Ci and C2 are positive constants depending only on q, ct', N and w^. Then, 
from (|8.84p and (|8.85p and by means of reflexivity, we get, passing to a subsequence 
if needed, that 

(8.89) pL^£° in 1^*^(0, T;A'^"(w')), as m ^ 00, 

(8.90) pL^e'' in L^+^(0,r;B''«(cj')), as m ^ 00. 

Now we pass to the limit m —> 00 in (|8.86p by using the convergence results (|4.40l) , 
(11311), and (|05)) . together with the identities u = |u|'^"2u ^^d G = u u, 

and also the convergence results (|8.89p and (|8.90p . Then we compare this limit 
equation with that one resulting from applying (|8.86p to the limit equation (|4.47p . 
considered for all (p S Cg°(a;^), with div (fi ~ 0, and with the aforementioned iden- 
tities. This procedure yields the existence of unique functions and p^ satisfying 
to (|8.84p - (|8.88p . Then, since p° and p^ are uniquely defined, we see that 

and 

Proceeding as in |24[ p. 126], letting a be such that 1 < a < 00, using the well- 
known local regularity theory, the compact imbedding W'^'"(cj) ^ W^''*(w) and 
Lebesgue's theorem of dominated convergence, we can prove that 

(8.91) P'L-^P'' strongly in L''{0,T;W^'''{uj)), as m ^ 00, 
where w is a fixed but arbitrary open bounded subset of il such that 

(8.92) UJ CC w' CC ri, with duj Lipschitz. 
Let us set now 

(8.93) v,„, :=u,„ +Vp,^^, 

(8.94) v:=u + Vp''. 

Then, combining (|8.86p with this same equation when we pass it to the limit m — > 
00, and using the definition of the distributive time derivative, we obtain 

(8.95) 

(v,„ - v)' = div (I Vu,„J''-2vu,„ - S) - a (|u,,„ r-2^,„^ _ |^|.-2^) 

inVlLUrp) 

-div(u,„ ® u,„,$,„(|u,„|) -u® u) - V(p^^ -/) 
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Proceeding as for (|5.51|) . attending to (|8.89p and observing that ro < r, we can 
prove that 

-v)' GL'■"(0,T;W-l'^°(w')) + L"'K)• 
Now we shall decompose the pressure term (p^ — p") in (|8.95[) into three new 
functions. For that, we need to invoke the following results, whose proofs follow 
immediately from \2A\ Lemmas 2.3 and 2.4]. 

Lemma 8.1. Let 1 < s < oo and fc G IN. 

(1) Then for every v* G ^Wq'* there exists a unique v £ Wq''^{U) such 
that 

f B"vD"ipdx= {v*,ip) V(^eCg°(a;), |a| = k. 

J u 

(2) In addition, if exists % G L*(ti;) such that 

■HD"(^dx V^gCS°(w), |a|=k. 



then 

||D%||l.(„) <C||H||l=m, 
where C is a positive constant depending on s and on the Calderdn-Zigmund 
inequality's constant. 

By a direct application of the first part of Lemma 18.11 attending to ()4.42p , (|4.43p 
and (|4.45p . and to the definitions of A'' {uj') and A3^^ (u)'), there exist unique 
functions 

(8.96) GL''(0,r;A''(w')), 

(8.97) pI^^ G L«^(0,r; A«^(w')), 

(8.98) GL^'(0,r;Wj''^'(^')), 
such that 

(8.99) / pl^A(j)dyidt^ I (jVue^l'^^Vue^ - S) : V^^dxdt, 

(8.100) / p^^^A(f>d:>idt = - (ue„ ® u,„$,„^(|u,,„|) - u«)u) : V^^dxdt, 

(8.101) / pl^A(l)dyLdt = a 1 (|ue„ |'""^Ue,„ - jul'^'^u) . V^dxrfi 

for all (j) G Cg°(a;^). In addition, by (|8.99p - (|8.10ip and a direct application of the 
second part of Lemma [5TT1 the following estimates hold: 

(8.102) IIpLIIl.'(c.^) <Cill|Vu,„J«-2yu,„-S||L,'(^^); 

(8.103) ||p2^J|^,«+2^^^^ <C2||u,„®u,„$,„(|u,„|)-u®u||^,«+2^^^^; 

(8.104) llVp3^||L.'(^g < C3|l|u,„r-2^,„ - |ur-2uilL.'(^g; 
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where Ci, C2 and C3 are positive constants depending on q' , q^^-^^ and a' , respec- 
tively, and on the Calderon-Zigmund inequahty's constant (C3 depends also on a). 
Next, testing (|8.95p by V0, with cj) £ C^{ujt), integrating over ujt and using p.l4p 
and ((T!M)) together with (fTM)) . and also the identities ((5?gg)) - (|8.10ip . we obtain 

Inserting this into (|8.95p . it follows that 
(8.105) 

(Ve,„ - V)' = div (|VU,„|9-2VU,„^ - S) - a (|Ue,„r-2ue„ - ju^-^u) 

- div(ue„ ig) Uc^$e„(|u,„J) - u(g) u) in V'iujT)- 

9. Definition of the irregularity regions 
Let us consider the following slight modification of the functions (j8.93p - (|8.94p 

(9.106) w,,„ (v,„ - v)xc.r = (u,,„ + Vp^^^ - (u + V/)) Xu^t, 

where x^t denotes the characteristic function of the set lot := uj x (0,T) and w 
satisfies to (|8.92p . Having in mind the extension of (|8.105p to R^+^, here we shall 
consider that 

(9.107) Ql^ := S - I Vu,„J«-2 vu,„ | + plj, 

(9.108) Ql^ :=u,„(8u,„$e„.(|ue„|)-u®u + p2^I, 

(9.109) :=a(|u,„r-2u,„ - \u\''~^u) +Vpl^, 

are extended from ujt to 1R,^+^ by zero. Now, since q < q* , we can use ()5.54p . with 
s = 7 = g, together with (|8.9ip . with a ^ q, to prove that 

(9.110) We,„^0 strongly in L'?(R^+i), as m ^ 00. 
Moreover, using (|4.34p and again (|8.9ip with a = g, we obtain 

(9.111) I|Vw,„J|l,(e«+i) < C. 
On the other hand, due to (|4.42p and (|8.102l) , we have 

(9.112) IIQLIIl.'(e«+i)<C. 

Moreover, ([STSSjl together with (|8.103p . and (j446l) together with (|8.104p . justify, 
respectively, that 

(9.113) Q2_^ ^ strongly in L'^^w^(R^+1), as m ^ 00, 

(9.114) qe,„^0 strongly in L'"'(R^+i), as m ^ 00. 

In order to define the irregularity regions of the admissible function that we shall 
test in (|8.105p . let us set 

(9.115) :=X*(|we„J), 

(9.116) .ge„. ■= M*{\S/^,J) + {M*i\QlJ))^ , 

(9.117) (Ai*(IQLl))^' 
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(9.118) :=(A^*(|q,,J))^, 

where M.* := M.t ° -A^x- Here Ait and Aix denote the Hardy-Littlewood maximal 
operators, which are defined, for some function / G LP(K,^+1) with 1 < p < oo, 
respectively by 

1 /■*+'^ 

7Wt(/)(x,i) sup — / |/(x,s)|ds, 

0<r<cx) ■ir Ji_j. 

MM)M--= sup -—^—j |/(y,,s)|dy, 

where Bu{x.) denotes the ball of centered at x and with radius R > 0. Then 
due to the boundedness of the operator A4* from L^'(R^+^) into L^'(K,^+^) for any 
p > 1 (see e.g. Stein [22l p. 5]), we obtain 

(9-119) \\feJ\l.<i(RN + r) < C||We„|jL<!(]Riv + i), 

(9.120) ||g.„J|L.(E«+i) < Ci||Vw,„||l,(r«+i) + C2\\QlJ0^^^^,y 

(9.121) ||X*(|Q? 1)11 "+2 <C||Q^ II iv+2 

(9.122) l|A^*(|q.„l)llL^'(R«+i) <^^l|q.„llL^'(R«+i)- 

Next, let be anyone of the functions inside the norms on the left-hand sides 
of (|9.119p - (|9.122p and let s be the respective Lebesgue exponent. Using (|9.11ip - 
(|9.114p . (|9.119p - (|9.122p and arguing as in [H P- 31], we obtain for j € IN 

IIj.™IIl=(R«+i) > 111(2) . inf t^Cn+i ({(x,t) G R^+i : b.„| > r}) . 



As a consequence, there exists Xm,j S 



22' 22 



such that 



(9.123) Cn+1 {(x,<) e R^+i : \j,J > X,n,j} < C2"^\~J^^ lb.„||LP(R«+i). 
Let us consider the following subsets of R^+i 

(9.124) F^,, := {(x,0 : |/.„J > A„,,}, 

(9.125) G™,, {(x,i) : |5,„| > A,„,,}, 

(9.126) H^^, := {(x,0 : \K_J > A™,,} ^ {(x,t) : X*(|QLI > A^^:)} , 

(9.127) {(x,i) : |*,„J > \^.,} ee {(x,t) : M*{\q.J > X^]} . 
Then, using (|9.123p in each case separately, we obtain 

(9.128) CN+iiF^^j)< G2-^A„'^.||/,„||l,(e«+i), 

(9.129) CN+i{G,nj)< G2"^A,7^-||.9,„||l,(k«-^i), 



(9.130) CN+i{Hm,j)< C2-^X„-,y" \\M*{\Cl:^An\,.N+2 



/ i\ N + 2 

(<?-l)9-2S- II A ^*n/^2 |\|| 



(9.131) £w+i(/™,,) < C2-^A-|f ||A^*(|q.„J)||L.'(K"+i)- 
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Now, since Xm,j G 



we observe that (|9.110|) . (|9.119p and (|9.128p on the 

one hand, (|9.113p . (|9.12ip and (|9.130p on the other, and yet (|9.114p . (|9.122p and 

(|9.13ip on another one, imply, respectively, 

(9.132) 

linisup£jv+i(^'mj) = 0, linisup£jv+i(-ffm.i) = 0, linisup£jv+i(-fmj) = 0. 

m^oo m— >-oo m— foo 

Moreover, since A^* is subadditive (see e.g. Stein jH]), we get from the definitions 
of G^j, Hm,j and Im,j in ([025| - (|9T27l) . using ((029l) - (|9l3T|l and ({91161) - (|9JT8| . 
that 

(9.133) G™j U U /™j D O, 
where 

O := {(x,t) € R^+i : M* (|Vw,,J) + (X* (|q1;2|) +X* (|q,,J)) < 4A„,,,} , 

(9.134) p™,, -.^ Af„7 
and, with the notations of (|9.107p - (|9.108p . 
(9-135) Qif :=Q,\„+Q^^. 
Setting 

(9.136) E,nj ■■= iF,n,j U G„j U U n wt, 

we can readily see that due to (|9.128p - (|9.13ip . (|9.119p - (|9.122p and to (|9.110p - 

dnm, 

(9.137) /:n+i(S„j) < oo- 
Moreover, due to (|9.133p . we have 

(9.138) {OLiU)nujT C E„i,j CUJT, 
where here U is the set F^.j defined in (|9.124p . 

10. Construction of a Lipschitz truncation 

We are now in conditions to define the truncation we shall consider here. Let us 
consider the following family of cubes 

(10.139) C,^--(x„,t„) {(y,s) e : dp^ ^ ((x„, t„), (y, s)) < r„} , 

where r„ > 0, n G IN and dp^ ■ is the metric defined by 



(10.140) dp^ ^ ((x„,i„), (y,s)) ;= max i |y - x„|, -y/p„^|^.|s - t„| 



By Theorem VI. 1.1])), there exists a Whitney covering of i?,„.j formed by the 
family of cubes (|10.139p - (|10.140p such that 

I J ^ 1 (^n ; ^n) — Ejji j . 

2 ' " 

neIN 

Moreover, by [2^ Section VI. 1.3], there exists a partition of unity n G M, 
associated to the Whitney covering (|10.139p - (|10.140p such that 

^Vfe = l m C,^',7-^(x„,t„)- 

feGM 
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We are now in conditions to define the Lipschitz truncation. Following |121 Section 
3] and [221 Chapter VI] , we define the Lipschitz truncation of w^^ subordinated to 
the Whitney covering (|10.139p - (|10.140p by 

We,„ in \ E„i,j 



(10.141) 7;„,(w.J:.^ E^^w. , ^ , ini.„. 



n=l 



The idea of this truncation, is to regularize the function by cutting off the re- 
gions Em,j of irregularity and then to extend this restricted function by the Whitney 
covering (|10.139p - p0.140p to the whole domain again. 

Now, let ^ G C'^{ujt) be a fixed cut-off function such that < ^ < 1 in wt and 
let us consider the following admissible test function for (|8.105p 

(10.142) :=e7;„ J (we„). 

In order to establish the main properties of the Lipschitz truncation (|10.14ip we 
are interested in, let 

(10.143) := supp^, ^ is the cut-off function of (|10.142p . 

Note that oj^ is strictly contained in lot, because < ^ < 1 in ut- Let also 
C^'^ . (w|n) be the space of all Lipschitz continuous functions with respect to the 
metric ()10.140p . From the definition of w^^ (see (|9.106p ). using (jOI]) and ([6?56l) 
together with ()8.9ip . with a = q, and (|8.88p . we can prove that 

We„ G L°°(0,r;L2(cj))nL«(0,r;Wi'«(w)). 

Then, owing to (|9.132p - (|9.138p . we can apply directly [HI Theorem 3.9, (i)-(iii)] to 
obtain: 

(10.144) TmA^,J^C"ti^i)^ 

with the norm depending on N, lu^, ||wj^J|li(_e,„ j), ||we^J|Li(Sr)j where CC 
u!t CC wt; 



(10.145) ||vr™,,(w„j|iL^(^£ ) < c (A™,, + p7n]jS;^;'d^eJ\i.HE^,,: 

(10.146) |17;„,,(w,„J||l^(^,) < C (l + p-|^.r^-2^||w,„J|L.(^„^ ^)) ; 

||(7;„j(w,„))' . (7;„,j(w,„J - (we,J)||Li(^enis ■) - 

(10.147) _^ _^ " 2 

CP,nj^N+l{E.mA (^A,„j + P™./^^^ ^_^£ ||We„J|Li(iS,„,,) j . 

In (|10.145p - p0.147p the constants C are distinct and depend only on N, and 

(10.148) ^p„,.4 dp^.ji^T^^r) > due to (|10.143p . 
Moreover, according to [121 Lemma 3.5] (see also [22l Section VI. 3]), 

(10.149) l|r™,,(w,„)||L.(^^) < C||w,,J|l.(^,) V s : 1 < s < oo, 
where C depends only on N. 
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11. Convergence of the approximated extra stress tensor 

Proceeding as for (|8.105|) . observing that now the functions are zero outside ujt 
and using the notations (|9.135p and (|9.106p - (|9.109p . we obtain 



(11.150) W^^ = -divQl'^ - c[^^ inX>'(wT) 



Here the distributive time derivative w' is such that 



(11.151) w^,^ e L''''(0,T;W-i^''«(a;)) (wt 



where tq is defined by ([5^. In fact, due to ([g^S]) . (|8.102p and (|8.103p on 

one hand, and due to (|4.46p and (|8.104p on the other, we can prove that 



(11.152) Ci\'^£-U'>{ujt) and qe„ G L""' (wt). 

As a consequence of (|11.152p i . diw^,]^ e L''"(0, T; W"i'''o(w)). Now, observing 
that, by virtue of (|10.144p - p0.146p and of the definition of ^, our admissible test 

function, defined in (|10.142p . <?!)™j e 17'^'{Q,T-MVY'° {uj)) n L'"(wt)- Then, from 
(|11.150p and pi.l5ip . we infer that 

(11.153) I (w^^ (<), 0,„,,(t)) dt= I Qi;2 : V(^,„,, d^dt - f • 0™,, dxdt. 

Jo J CUT J '^T 



On the other hand, owing to (|9.132p - (|9.138p and, in addition, to (|11.152p - pi.l53p . 
we can apply [121 Theorem 3.9, (iv)] to prove that for every ^ e Cjf (w^) 



T 

(11-154) (ir™,j(w,„)i'-2w,,„-7;,j(we„))e'dxdi- 

(7^ij(We„))' • {%n,j{^e„J - ^^e^)^dxdt. 



Note that the proof of (|11.154p is done in [TH p. 23] for q^^ ee in (|11.153p . But 
taking into account (|11.15ip . the proof of [I^ Theorem 3.9, (iv)] can be repeated 
almost word by word in our case. 
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Now, gathering ()11.153p and (|11.154|1 . and expanding the notations ()9.135|l and 
([9l06| - ((9T^ . we obtain 

f (|Vu,,J«-2vu,„ - S) : V(r™,j(w,,J)ec;xdi = 

J LOT 

+ f (S - |Vu,„|9~2y„^^^ . r,n,j{w,J ® dxdi 

J LOT 

+ / (Ue„ «) Ue„$e„(|Ue„|) - U(g) U) : V (7;„j(We,„)0 rfxdt 

J LOT 

J LOT 
J LOT 

^^^■^^^^ + f pi^div(r™,,(w,,j)edxdt 

J LOT 

+ / pl^AYv{Tm,j{'We„,)C)dyidt 

J LOT 

' / ^vl-^-r^A^cJidy^dt 

J LOT 

+ / T'nji^im) ■ - T^n.ilWe^)) £, dxdt 

J E,„ , 



Ji + J2 + J?, + Ji + J5 + ■h + -h + J& + Jg- 



We claim that, for a fixed j, 



(11.156) limsup 



/" (iVu.^J'-^Vu,,^ - S) : V(T^,,(w,„))erfxdt 

^ LOT 



< C2- 



To prove this, we will carry out the passage to the limit m — > 00 in all absolute 
values \ Ji\, i — 1, . . . , 9. 

• limsup„^^(|Ji| + IJ4I) = 0. Due to ((i:!^ and ^JU^, S- | Vu^^ I'-^Vu.^ 
and pI are uniformly bounded in L'^ (wt)- Then, using Holder's inequality and 
P0.149P together with (|9.106p . led us to 

I Jil+I J4I < Ci \\%nA^.JL.i^r) ^ ^2 (||u,,„ - u||l,(^,) + ||v(pt -?5")IIl^(c..)) ■ 

The assertion follows by the application of (|5.54p with s = 7 = q and (|8.9ip with 
a ~ q, and observing that always q < q* for any q > I. 

• limsup,„_j.o2(| J2I + l^el) = 0. In fact, by Holder's inequality, 

I J2I + I Jel < |lUe„ ® U,„$,„,(|U,,J) - U® u||l1(^^)||V (7;„jK„)C) + 

lbLllLi(^T)l|div (7mj(We,„)0 IIl-(c^|,)- 

Then, using Holder's inequality again and due to (|5.55p and (|8.103p . we get 
limsup^_^gc(| J2I + I Jgl) = if both second multiplying terms on the right-hand 
side of the above inequality are finite. Indeed, by the application of (|10.145p and 
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()10.146p together with (|9.106p . we get 



Ci A, 



C2 1 



V T.l 



From (|4.40p and H8.9ip . the last with a = q, v^^ — v is uniformly bounded in 
L^(_Emj). On the other hand, for a fixed j e IN, the sequence X^.j lies in the 



interval 



and, as a consequence, the sequence p„i,j = A^j J is uniformly 



bounded from above, which, in turn, by (|10.148p and (|10.140p . implies 

inf 6^ . ^5 > 0. 

Analogously, we prove that also ||div (7^j(we^) ^) ||loo(^« ) is finite. 
• limsup„j_j.o2(| J3I + IJ7I) = 0. By Holder's inequality and (|8.104p 



|J3| + |J7|< a|||u. 



'ue„ - |u|" 'ui|Li(^^) + ||Vp3^J|Li(„^)) I|7;„,j(w,„)||l^(^5^) 



Arguing as in the previous case, we can show that, for each j e IN, Tm,j{^e„-,) is 
uniformly bounded in L°°(cjy). Then, by the application of (|4.46p . it follows that 
limsup,„^^(| J3I + I J7I) = 0. 

• limsupjjj^g^d J5I + 1 J9I) < C2~~i . By the definition of the Lipschitz truncation 
(see (|10.14ip ) together with the fact that divw^^ = (see (|9.106p ). we can write 



J5 = / p\^d\YTm,]{y^t^)d-iidt. 



Next we use Holder's inequality, and (|8.102p together with (|10.145p . Then, 

arguing as in the cases for | Ji| + | J4I and | J2I + | Jel, we have 



limsuplJsl < Cilimsup||V7^„j(we^)||L,. £p,£;^^ .) 

< Ci limsup£Ar+i(i;™j)9 ||V7^j(vife„J||Loo(^« ) 



< C2 lim sup 



< C2 lim sup ( Cn+1 [Em.,] ) " A„ 



Next, by the definition of E^^^ (see (|9.136p ) and using fJ.132p . (|9.129p . (|9.120p and 
(|9.11ip - (|9.1f2p by this order, we get limsup^^^ l-M < C2~i . 
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For Jg, we have by using (|10.147p together with the definition of pmj (see (|9.134p ) 
and arguing as we did above for | Jsj, 



hm sup I Jg I < C Hm sup 



|ve„ - v||Li(i-;,„.,) 



< CHmsup ( )£Ar+i(£;,„,j ) 



m— fcxD 



Then observing that g > 1, it follows that limsup„ 



iJgl <C2- 



Throughout the above bullets, we have proven the claim (|11.156p is true. On the 
other hand, arguing as we did for | J5I, we can prove also that, for a fixed j. 



(11.157) limsup 



[ (|Vu,,„ r-2Vu,„^ - S) : V7;„,,(w,,J fdxdi 

Je„, i 



In consequence, from the definition of Tmj (see (|10.14ip ). ()11.156p and (|11.157p . 
we prove that 



(11.158) limsup 



/ (I Vu,„ 1''"^ Vu,,„ - S) : Ww,^ ^ dxdt 

J UJT\Em,j 



Using the definition of w^^ (see (|9.106p i and the strong convergence property of 
(see (|8JT|) ). it can be derived from (|11.158p that 



(11.159) limsup 



f (I Vu,,„ 1*^"^ Vu,„ - S) : V(u,„ - u) ^ dxdt 

J LJT\Em.j 



< C2- 



Now, by (|11.159p and (|9.132p . for each j S IN we can find a number rrij € IN such 
that 



JuiT\Emj ,j 

(11.161) CN+i{F,n„j) < C2-\ 

(11.162) /:w+i(Jf,n,j) <C2-^ 

(11.163) /:7v+i(/m,,j) < C2-J. 

Setting := ^Xajr\s,7i j > where Xtjy\£;^ ■ denotes the characteristic function of the 
set ujT\Em,,j, it can be proved (c/. [H PP- 36-37]), using (pi^ . pi.l6ip - (|11.163p 
and the fact that Xm^^j > 1, that 

(11.164) — > ^ a.e. in lut as j — !• od. 

From pi.l64p . (|4.40p and (|4.42p . we have, by appealing to Lebesgue's theorem of 
dominated convergence, that 

(11.165) Vu^j ^ Vu^ strongly in 'L'^{ll!t), as j — oo, 

(11.166) S^j^S^ strongly in L?'(a;T), as j ^ 00. 



EXISTENCE FOR THE GENERALIZED NAVIER-STOKES EQUATIONS WITH DAMPING 23 



Then, from (|11.160p and ()11.165p - pi.l66p . and appealing once more to dMQ]), 
()4.42p and Lebesgue's theorem of dominated convergence, we obtain 

(11.167) lim / |Vue„J«^j dxrff= / S : \/u£_ dxdt. 

Finally, taking into account ^M^, (|11.164p and (|11.167p . we can apply the 

local Minty trick (c/. [H Lemma A. 2]) to establish that = jVul'-^Vu^ a.e. 
in UJt- Due to the arbitrariness of ^, S = |Vu|'^^^Vu a.e. in ljt and the proof of 
Theorem 12. H is concluded. ■ 

12. Remarks 

In Theorem 12.11 we have proved the existence of weak solutions, in the sense of 
Definition \TT\ to the problem prTjl - PTi)) for any 

It is only left open the case oil < q < ior N > 3. But with the methods at our 
present disposal it seems to be very difficult to prove this case, because the compact 
imbedding W^''(il) L^(il), which holds only for q > 7^^, is fundamental in 

many steps of our proof. 

The result established in Theorem 12.11 is still valid if we consider an extra stress 
tensor with a g-structure satisfying to general growth and coercivity conditions. 
Indeed the proof still holds with minor changes if we assume that the diffusion term 
I Vu|'^~^Vu in (|1.2p is replaced by a tensor T = T(x, t, D) (D is the symmetric part 
of Vu) satisfying to 

• |T(x,t, A)| < Ci|A|9-i + /i for ah A in M^'y„, for a.a. (x,t) in Qt and 

for any function /i in L'? (Qt) with /i > 0; 

• T(x,t,A) : A > C2|A|« - /z for all A in M^y,„, for a.a. (x, i) in Qt and 
for any function /2 in L^{Qt) with /2 > 0; 

where Ci and C2 denote positive constants and M"yjj-j is the vector space of all 
symmetric n x n matrices, which is equipped with the scalar product A : B and 
norm |A| = vATA. 

It is possible to consider unbounded domains with no restriction on the size 
and shape of Jl. In this case, proceeding as in Section 3], we can prove the 
regularized problem (|3.14p - (|3.17p has a unique weak solution for such ft. As a 
consequence the original problem has a solution for these domains as well. 

The uniqueness of weak solutions is, as is well known, an open problem for 
the generalized Navier-Stokes problem (without damping) for values oi q < 2. 
By adapting [T71 Theoreme 2.5.2], we can prove the weak solution to the prob- 
lem (|l.ip - (|1.4p is unique under more restrictive conditions that we have needed to 
prove the existence. In fact, assuming that q > replacing the diffusion term 

div (|Vu|^^^Vu) in p.2p by div (|Vu|'~^Vu) + Au and having in mind the damp- 
ing term satisfies to (j2.5p , it is possible to prove the uniqueness of weak solution to 
this modified problem pTT|) - ([0)) (see [TBI Theorem 2]). 

A completely different issue, is the important question about the qualitative 
properties of the weak solutions to the problem (|l.ip - (|1.4p . In this scope we are 
mainly interested in the extinction in a finite time, once that the confinement of 
the weak solutions in a space domain is a much more delicate matter and remains 
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an open problem, with the exceptions of the stationary Stokes and Navier-Stokes 
problems (c/. [l]-[4] and [E]). Proceeding as in [HI Theorem 3], letting u be a 
weak solution to the problem (|l.ip - (jl.4p in the sense of Definition l2.1l and assuming 
that (|2.10p and one of the following conditions hold: 

(1) q<2- or 

(2) 1< cr < 2; 

then we can prove the following assertions are true with minor changes in the proofs: 

• If f = a.e. in Qt, then there exists i^^^ > such that u(x,t) = a.e. in 

and for all t > t*^,.y for i = 1, 2; 

• Let f 7^ and assume that exist positive constants ej^) and (positive) times 
ij'^j, for i = 1, 2, such that, for almost all t G [0,T], 



or 




Then there, there exist positive constants e^^ such that u = a.e. in il 

and for all t > t^^^ provided < ej-.j) < e^-y for i = 1, 2. 

Note that the subscripts are used above in the sense to relate each result to the 
different condition (i) written before. 
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